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Recently, it was observed that water droplets suspended in a nematic liquid crystal form linear 
chains (Poulin et al, Science 275, 1770 (1997)). The chaining occurs, e.g., in a large nematic drop 
with homeotropic boundary conditions at all the surfaces. Between each pair of water droplets a 
point defect in the liquid crystalline order was found in accordance with topological constraints. 
This point defect causes a repulsion between the water droplets. 

In our numerical investigation we limit ourselves to a chain of two droplets. For such a complex 
geometry we use the method of finite elements to minimize the Frank free energy. We confirm an 
experimental observation that the distance d of the point defect from the surface of a water droplet 
scales with the radius r of the droplet like d » 0.3 r. When the water droplets are moved apart, 
we find that the point defect does not stay in the middle between the droplets, but rather forms a 
dipole with one of them. This confirms a theoretical model for the chaining. Analogies to a second 
order phase transition are drawn. We also find the dipole when one water droplet is suspended in 
a bipolar nematic drop with two boojums, i.e., surface defects at the outer boundary. Finally, we 
present a configuration where two droplets repel each other without a defect between them. 

PACS numbers: 77.84.Nh, 61.30.Cz, 61.30.Jf 



I. INTRODUCTION 

Liquid crystal research is, without any doubt, strongly 
driven by its importance for technological applications, 
namely display devices. On the other hand, topologi- 
cal point and line defects in the liquid crystalline order, 
which occur as a consequence of the broken rotational 
symmetry of the isotropic space jl|-[5| , attract a lot of at- 
tention since they are fascinating by themselves. They, 
furthermore, influence the static and dynamic properties 
of the system under investigation. Conventional displays, 
like twisted nematic or surface-stabilized ferroelectric 
liquid crystal cells J?]], possess a simple geometry, where 
the liquid crystal is enclosed between two parallel glass 
plates. With the introduction of polymer-dispersed liq- 
uid crystals as electrically controllable light shutters , 
an extensive study of liquid crystals confined to complex 
geometries, like drops in a polymer matrix or a random 
porous network in silica aerogel, was initiated |@. 

This article presents a numerical investigation of the 
interesting inverse problem, which is posed by particles 
suspended in a nematic solvent. In a nematic liquid crys- 
tal the molecules align on average along a unit vector n, 
called director. The energetic ground state is a uniform 
director field throughout space. Due to the anchoring of 
the molecules on the surface of a particle the surrounding 
director field is elastically distorted. In addition to con- 
ventional van der Waals forces, screened Coulomb and 
steric interactions JllJ] , the elastic deformation of the di- 
rector field mediates a further interaction between the 



particles. Point defects give rise to a short-range repul- 
sion. 

Already in 1970, Brochard and de Gennes studied a 
suspension of magnetic grains in a nematic liquid crys- 
tal and determined the director field far away from a 
particle fll2|| . A bistable liquid crystal display was intro- 
duced based on a dispersion of agglomerations of silica 
spheres in a nematic host Chains and clusters 

were observed in the dispersion of latex particles in a 
lyotropic liquid crystal 16 The radii of the parti- 
cles were 60 and 120 nm. Therefore, details of the di- 
rector field could not be resolved under the polarizing 
microscope. Terentjev et al. investigated the director 
field around a sphere by both analytical and numerical 
methods [|l9|,^0). For homeotropic boundary conditions 
of the director they considered the Saturn-ring configu- 
ration where the sphere is surrounded by a —1/2 discli- 
nation ring. The long-range quadrupolar interaction of 
such objects was determined by Ramaswamy et al. |2l]] 
and Ruhwandl and Terentjev p2| . Even two parallel 
plates, immersed into the isotropic phase close to the 
isotropic-nematic phase transition, interact as a result of 
a surface- induced nematic order |23|. 



A. Inverted Nematic Emulsions 

The present article addresses recent work on in- 
verted nematic emulsions where surfactant-coated water 
droplets are dispersed in a nematic solvent 24 2^] . The 
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advantage of such a system is that the particles are easily 
observable by polarizing microscopy since the size of the 
water droplets is of the order of a micron. Furthermore, 
the anchoring of the liquid crystal molecules on the sur- 
faces of the droplets are controllable via the surfactant. 
Understanding the suspension of such large particles will, 
without any doubt, help to clarify the observations of dis- 
persions of smaller objects. The most striking feature, 
from a theoretical point of view, is that inverted emul- 
sions provide an ideal laboratory for the investigation of 
topological defects. In inverted emulsions point defects, 
also called hedgehogs, occur. They carry a topological 
charge q specifying the number of times the unit sphere 
is wrapped by the director n on any surface enclosing 
the defect core. Since all properties of the nematic phase 
are invariant under the inversion n — > — n, q is always 
positive, and two hedgehogs with respective charges qi 
and (?2 "add up" to a total topological charge of either 
5l + 92 or fa - q 2 \ f§. 

When a water droplet with homeotropic, i.e., perpen- 
dicular anchoring of the director at its surface is placed 
into a uniformly aligned nematic liquid crystal, a hyper- 
bolic hedgehog is nucleated (see Fig. |l|) . The droplet and 
its tightly bound companion defect provide a key unit 
to understand inverted nematic emulsions p^j2^ , |25| , p7t . 
We will call it a dipole for short because of its dipolar 
symmetry. Both the water droplet and the hyperbolic 
hedgehog carry a topological charge +1 which "add up" 
to the total charge of the dipole. A phenomenological 
theory predicts that the dipole couples to a strong splay 
deformation in the director field plp^ ]. Two parallel 
dipoles are attracted via a long-range dipolar interaction 
which explains the observed chaining of water droplets 
P^ , p5|j27| . A quantitative confirmation of the dipolar in- 
teraction was presented recently |28[ ] . 

In our numerical investigation we consider multiple ne- 
matic emulsions where the nematic host fluid containing 
the water droplets is confined to large nematic drops, 
which, in turn, are surrounded by water. The water 






FIG. 1. The water droplet and its companion hyperbolic 
hedgehog form a dipole. Both the droplet and the hedgehog 
carry a topological charge +1, and the total charge of the 
dipole is = |1 — 1|. 



FIG. 2. Scenario to explain the chaining of water droplets 
in a large nematic drop. The right water droplet and its com- 
panion hyperbolic hedgehog form a dipole, which is attracted 
by the strong splay deformation around the droplet in the 
center (left picture). The dipole moves towards the center 
until at short distances the repulsion mediated by the point 
defect sets in (right picture). 

droplets also form chains with a hyperbolic he dgehog sit- 
uated between two droplets (see Fig. ^, right) 24 2j|. In 
the experiment it was found that the distance d of the 
point defect from the surface of a water droplet scales 
with the radius r of the droplet like d ss 0.3 r In 
the following we will call this relation the scaling law. A 
possible explanation for the chaining, as in a uniformly 
aligned sample, could be the presence of a dipole. One 
water droplet fits perfectly into the center of a large ne- 
matic drop, which has a total topological charge +1. Any 
additional water droplet has to be accompanied by a hy- 
perbolic hedgehog in order not to change the total charge. 
If the dipole forms (see Fig. |^, left), it is attracted by the 
strong splay deformation in the center |27j , until a short- 
range repulsion mediated by the defect sets in p4,p5[ . 



B. Main Results 

By our numerical investigations we confirm the scaling 
law, for which we find d = (0.325 ± 0.025) r, and discuss 
the influence of the outer boundary of the large drop. We, 
furthermore, show that in nematic drops the dipole also 
exists, confirming the scenario of Fig. |^. When the two 
water droplets in the right picture of Fig. || are moved 
apart symmetrically about the center of the large drop, 
the dipole forms via a second order phase transition. We 
also identify the dipole in a bipolar configuration which 
occurs for planar boundary conditions at the outer sur- 
face of the nematic drop. Two boojums, i.e., surface 
defects appear [p9|-[3lf , and the dipole is attracted by the 
strong splay deformation in the vicinty of one of them 
pi] , p5 27 1 . Besides the dipole we find another stable con- 
figuration in this geometry, where the hyperbolic hedge- 
hog sits close to one of the boojums. Finally, we show 
that water droplets can repel each other without a hy- 
perbolic defect placed between them. 

To minimize the Oseen-Zocher-Frank free energy p3| l 
for different positions of the water droplets we employ the 
method of finite elements 32 1 , which is most suitable for 
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FIG. 3. Geometry parameters for two water droplets with 
respective radii n and in a large nematic drop with radius 
rz- The system is axially symmetric about the z axis, and 
cylindrical coordinates p, z are used. The coordinates zi, z%, 
and Zd are the respective positions of the two droplets and 
the hyperbolic hedgehog. The two distances of the hedgehog 
from the surfaces of the droplets are di and d^. 

non-trivial geometries. Since we are not aware that this 
method has been used before in liquid crystal research, 
we provide a more detailed description of it in Sec. || and 
in the Appendix. We hope to stimulate further applica- 
tions of the method of finite elements to liquid crystals 
in complex geometries. Before we discuss our results in 
Sec. [II we define the geometries and present numerical 
details in Sec. |fi|. 



II. GEOMETRIES AND NUMERICAL METHODS 

In our article we numerically investigate two particular 
geometries of axial symmetry. The first problem is de- 
fined in Fig. [| We consider two spherical water droplets 
with respective radii r% and T2 in a large nematic drop 
with radius r 3 . The whole system possesses axial symme- 
try, so that the water droplets and the hyperbolic hedge- 
hog, indicated by a cross, are located always on the z 
axis. We employ a cylindrical coordinate system. The 
coordinates z-y, z%, and z,i denote, respectively, the posi- 
tions of the centers of the droplets and of the hyperbolic 
hedgehog on the z axis. The distances of the hedgehog 
from the surfaces of the two water droplets are, respec- 
tively, d\ and di- Then, the quantity d\ + di means the 
distance of the two small spheres, and the point defect 
is situated in the middle between them if d\ = d%. We, 
furthermore, restrict the nematic director to the (p, z) 
plane, which means that we do not allow for twist defor- 
mations. The director is expressed in the local coordinate 



basis (e p ,e z ), 

n(p, z) = sin B(p, z)e p + cos e(p, z)e z 



(1) 



where we introduced the tilt angle B. It is always re- 
stricted to the range [— 7r/2,7r/2] to ensure the n — > — n 
symmetry of the nematic phase. At all the boundaries 
we assume a rigid homeotropic anchoring of the direc- 
tor, which allows us to omit any surface term in the free 
energy. 

In the second problem we have only one water droplet 
inside a large nematic drop. We use the same coordinates 
and lengths as described in Fig. but omit the second 
droplet. The anchoring of the director at the outer sur- 
face of the large nematic sphere is rigid planar. At the 
surface of the small sphere we again choose a homeotropic 
boundary condition. 

We write the Oseen-Zocher-Frank free energy density 
Q for the director field ([[]) in units of K^/a 2 , where 
-K33 is the bend elastic constant and a the characteristic 
length scale of our system, typically several microns: 



/(e) = 
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(2) 



All lengths are given relative to a, and K\\ = K11/K33 is 
the scaled splay constant. The indices p and z indicate, 
respectively, partial derivatives d/dp and d/dz. 

In order to find the equilibrium director field the re- 
duced free energy in units of -K33C1, 



F 



d 3 xf(Q) 



(3) 



is minimized on a grid. Because of the nontrivial geom- 
etry of our problem, we decided to employ the method 
of finite elements |3^] , where the integration area is cov- 
ered with triangles. We note, that for simpler geometries 
the method of finite differences is usually used, where the 
grid is defined by the coordinate lines. We construct a 
net of triangles by covering our integration area with a 
hexagonal lattice with lattice constant b. Vertices of tri- 
angles that only partially belong to the integration area 
are moved onto the boundary along the radial direction 
of the appropriate sphere. As a result, extremely obtuse 
triangles occur close to the boundary. We use a relax- 
ation mechanism to smooth out these irregularities. The 
final triangulation is shown in Fig. || In the area be- 
tween the small spheres, where the hyperbolic hedgehog 
is situated, the grid is further subdivided to account for 
the strong director deformations close to the point defect. 

The reduced free energy F in units of is written 
as a function of the tilt angles Qi (i = 1, . . . , n) at the n 
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FIG. 4. Triangulation of the integration area (lattice con- 
stant: b — 0.495). Between the small spheres a refined net of 
triangles is chosen. 

vertices of the grid: F = F[®], where = (9i, . . . , 9„). 
For further explanations how F[&] is constructed within 
the method of finite elements we refer the reader to the 
Appendix. To find a minimum of the free energy we start 
with a configuration that already possesses the hyper- 
bolic point defect at Zd and let it relax via the standard 
Newton-Gauss-Seidel method |341: 



e 



old 



dF[&] ld 2 F[&] 



del 



(4) 



The derivatives F with respect to Qi are performed nu- 
merically. Although the point defect is not fixed to its 
original position during the relaxation it never moves into 
the global minimum, but rather stays where we place it. 
We use this fact to plot energy curves versus the position 
of the point defect. Other starting configurations, e.g., a 
uniform director field throughout the integration area or 
randomly oriented directors on the grid, do not produce 
the hyperbolic hedgehog, but instead relax into a config- 
uration with a higher energy, which we briefly investigate 
in subsection III E . 



Integrating the free energy density over one triangle 
yields a line energy, i.e., an energy per unit length. As 
a rough estimate for its upper limit we introduce the 
line tension /; = (Kn + l)/2 of the isotropic core of a 
disclination pq |. Whenever the numerically calculated 
local line energy is larger than // , we replace it by fi . As 
a result, the hyperbolic point defect is stabilized against 
opening up to a disclination ring. 

All our calculations are preformed for the nematic liq- 
uid crystal pentylcyanobiphenyl (5CB), for which the ex- 
periments were done p^j25[ |. It has a bend elastic con- 
stant -K33 = 0.53 x 10 dyn and a scaled splay elastic 



constant Kn = Kn/K^ = 0.79. The experimental ra- 
tio r^/r 1/2 of the radii of the large and small drops is in 
the range 10-50 ||§|2|]. The difficulty is that we want to 
investigate details of the director field close to the small 
spheres which requires a fine triangulation on the length 
scale given by ri/2. To keep the computing time to a 
reasonable value we choose the following lengths: r% = 7, 
r i/2 = 0-5 ... 2, and b = 0.195 for the lattice constant of 
the grid. In addition, we normally use one step of grid 
refinement between the small spheres (geometry 1) or be- 
tween the small sphere and the south pole of the large 
nematic drop (geometry 2). With such parameters we 
obtain a lattice with 2200-2500 vertices. 



III. RESULTS AND DISCUSSION 

In this section we discuss the results from our numer- 
ical investigation. We will mainly address the following 
questions: 

1. Can we confirm the scaling law d]y 2 ~ ^■^ r i/2i which 
was observed in experiment, by varying the different 
lengths in our geometry? 

2. Does the dipole formed by a water droplet and its com- 
panion hyperbolic defect, which is well-established for a 
director field uniformly aligned at infinity, also have a 
meaning in more complex geometries? 

3. Is the hyperbolic hedgehog necessary to mediate a 
repulsion between the water droplets? To answer this 
question we will investigate a configuration with discli- 
nation rings around the small spheres. 



A. Scaling Law 

We start with the first question. In Fig. || we plot 
the reduced free energy F as a function of the distance 
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FIG. 5. The free energy F as a function of the distance 
di + c?2 between the small spheres which are placed symmet- 
rically about z = (n = ri = 1). Curve 1: Zd = 0, curve 2: 
position Zd of the defect can relax along the z axis. 
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d\ + d 2 between the surfaces of the small spheres, which 
are placed symmetrically about the center, i.e., z 2 = —z±. 
Their radii are r\ = r 2 = 1- Curve 1 shows a clear 
minimum at di+d 2 ~ 0.7, the defect stays in the middle 
between the two spheres at z^ = 0. In curve 2 we move 
the defect along the z axis and plot the minimum free 
energy for each fixed distance d\ + d 2 • It is obvious that 
beyond d\ + d 2 = 2 the defect moves to one of the small 
spheres. We will investigate this result in more detail in 
the following subsection. 

In Fig. U we take three different radii for the small 
spheres, r\ — r 2 — 0.5,1,2, and plot the free energy 
versus d\jr\ close to the minimum. Recall that d\ is 
the distance of the hedgehog from the surface of sphere 
1. Since for such small distances d\ + d2 the defect al- 
ways stays at Zd = 0, i.e., in the middle between the 
two spheres, we have d\jr\ — d 2 /r 2 . The quantity -F m i n 
refers to the minimum free energy of each curve. For each 
of the three radii we obtain an energetically preferred dis- 
tance d\/r\ in the range of [0.3,0.35], which agrees well 
with the experimental value of 0.3. Why does a scaling 
law of the form di/ 2 = (0.325 ± 0.025) r 1 / 2 occur? If the 
small spheres are far away from the surface of the large 
nematic drop, its finite radius r% should hardly influence 
the distances d\ and d 2 . Then, the only length scale in 
the system is r\ = r 2l and we expect d\j 2 cx r 1 / 2 . How- 
ever, in Fig. |^ the influence from the boundary of the 
large sphere is already visible. Let us take curve 2 for 
spheres with radii ri/ 2 — 1 as a reference. It looks pretty 
symmetric around d\jr\ = 0.35. The slope of the right 
part of curve 3, which corresponds to larger spheres of 
radii ri/ 2 = 2, is steeper than in curve 2. Also, the loca- 
tion of the minimum clearly tends to values smaller than 
0.3. We conclude that the small spheres are already so 
large that they are strongly repelled by the boundary of 
the nematic drop. On the other hand, the slope of the 
right part of curve 1, which was calculated for spheres of 
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FIG. 6. The free energy F — F m \ n as a function of 
di/ri = d2/r2- The small spheres are placed symmetrically 
about z — 0. Curve 1: ri = r% = 0.5, curve 2: ri = r% = 1, 
and curve 3: n = r-z — 2. 



radii ri/ 2 — 0.5, is less steep than in curve 2. This leads 
to the conclusion that the boundary of the nematic drop 
has only a minor influence on such small spheres. 

When we move the two spheres with radii r 1 / 2 = 1 
together in the same direction along the z axis, the de- 
fect always stays in the middle between the droplets and 
obeys the scaling law. We have tested its validity within 
the range [0,3] for the defect position Zd- Of course, the 
absolut minimum of the free energy occurs in the sym- 
metric position of the two droplets, z 2 = —z\. 

We further check the scaling law for n ^ r 2 . We inves- 
tigate two cases. When we choose r\ — 2 and r 2 = 0.6, 
we obtain di/ 2 ~ 0.3 r^- In the second case, r\ = 2 
and r 2 = 1, we find d\ ~ 0.37 r± and d 2 « 0.3 r 2 . As 
observed in the experiment the defect sits always closer 
to the smaller sphere. There is no strong deviation from 
the scaling law d Y j 2 — (0.325 ± 0.025) r\i 2 , although we 
would allow for it, since n ^ r 2 . 

B. Identification of the Dipole 

The second question is if the dipole formed by one 
water droplet and a companion hyperbolic point defect 
has a meaning in our geometry. To answer this question 
we place sphere 2 with radius r 2 — 1 in the center of the 
nematic drop at z 2 = 0. Then, we determine the energet- 
ically preferred position of the point defect for different 
locations Z\ of sphere 1 (r\ — 1). The position of the 
hedgehog is indicated by A = (d 2 — d\)/{d\ + d 2 ). If the 
defect is located in the middle between the two spheres, 
A is zero since d\ — d 2 . On the other hand, if it sits at 
the surface of sphere 1, A is one since di = 0. In Fig. 
we plot the free energy F versus A. In curve 1, where the 
small spheres are farthest apart from each other {z\ = 5), 
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FIG. 7. The free energy F as a function of 
A = (d,2 — di)/(d\ + g^)- Sphere 2 is placed at Z2 = 0. The 
position z\ of sphere 1 is the parameter. Curve 1: z\ — 5, 
curve 2: z\ = 4, curve 3: z\ = 3.5, and curve 4: zi = 3. The 
radii are n = ri = 1. 
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we clearly find the defect close to sphere 1 . This verifies 
that the dipole is existing. It is thermally stable, since a 
rough estimate of the thermally induced mean displace- 
ment of the defect yields 0.01 |27]]. When sphere 1 is 
approaching the center (curve 2: z% = 4 and curve 3: 
Z\ = 3.5), the defect moves away from the droplets until 
it nearly reaches the middle between both spheres (curve 
4: z\ = 3). This means, the dipole vanishes gradually 
until the hyperbolic hedgehog is shared by both water 
droplets. 

An interesting situation occurs when sphere 1 and 2 
are placed symmetrically about z = 0. Then, the defect 
has two equivalent positions on the positive and negative 
part of the z axis. In Fig. |^ we plot again the free 
energy F versus the position A of the defect. ^From 
curve 1 to 3 (zi = zi = 4, 3, 2.5) the minimum in F 
becomes broader and more shallow. The defect moves 
closer towards the center until at around z\ = —z% = 2.3 
(curve 4) it reaches A = 0. This is reminiscent to a 
symmetry-breaking second order phase transition P]35| 
which occurs when, in the course of moving the water 
droplets apart, the dipole starts to form. We take A as 
an order parameter, where A = and A / describe, 
respectively, the high- and the low-symmetry phase. A 
Landau expansion of the free energy yields 



(1) 



(2) 



F(A) = F (zi) + a [2.3 - z x \A 2 + c( Zl )A 4 



(5) 



where z\ = —Z2 plays the role of the temperature. Odd 
powers in A are not allowed because of the required sym- 
metry F(A) — F(—A). This free energy qualitatively de- 
scribes the curves in Fig. ||. It should be possible to ob- 
serve such a second order phase transition with a method 
introduced recently by Poulin et al. |28| to measure dipo- 
lar forces in inverted nematic emulsion. When the two 
small droplets are filled with a magnetorheological fluid 




A = (d 2 -d 1 )/ {d x + d 2 ) 



FIG. 8. The free energy F as a function of 
A = (d,2 — di)/(di + (fe). The small spheres are placed sym- 
metrically about 2 = 0. Curve 1: z\ = —22 = 4, curve 
2: zi = —Z2 — 3, curve 3: z\ = — 22 = 2.5, curve 4: 
z\ — —Z2 ~ 2.3, curve 5: z\ — —22 ~ 2. The radii are 
ri = r 2 = 1. 




FIG. 9. Planar boundary conditions on the outer surface 
of the large sphere create boojums, i.e., surface defects at the 
north and the south pole. A water droplet with homeotropic 
boundary conditions nucleates a hyperbolic hedgehog. Two 
configurations exist that are either stable or metastable de- 
pending on the position of the water droplet; (1) the dipole, 
(2) the hyperbolic hedgehog sitting at the surface. 

instead of pure water, a small magnetic field of about 
100 G, applied perpendicular to the z axis, induces two 
parallel magnetic dipoles. Since they repel each other, 
the two droplets are forced apart. When the magnetic 
field is switched off, the two droplets move towards each 
other to reach the equilibrium distance. In the course of 
this process the phase transition for the dipole should be 
observable. 



C. The Dipole in a Bipolar Configuration 

It is possible to change the anchoring of the direc- 
tor on the outer surface of the large nematic drop from 
homeotropic to planar by adding some amount of glycerol 
to the surrounding water phase |24|. Then, the bipolar 
configuration for the director field appears 1 30 3l[] , where 
two boojums pS[|, i.e., surface defects of charge 1 are situ- 



ated at the north and south pole of the large nematic drop 
(see configuration (1) in Fig. ^|). The topological point 
charge of the interior of the nematic drop is zero, and ev- 
ery small water droplet with homeotropic boundary con- 
dition has to be accompanied by a hyperbolic hedgehog. 
In the experiment the hedgehog sits close to the water 
droplet, i.e., the dipole exists and it is attracted by the 
strong splay deformation close to the south pole |24(] , as 
predicted by the phenomenological theory ]24|j27j ] . 

A numerical analysis of the free energy F is in agree- 
ment with experimental observations but also reveals 
some interesting details which have to be confirmed. In 
Fig. 10 we plot F as a function of the position z\ of the 
small water droplet with radius r± = 1 . The diagram con- 
sists of curves (1) and (2), which correspond, respectively, 
to configurations (1) and (2) in Fig. The free energy 
possesses a minimum at around z\ = —5.7. The direc- 
tor field assumes configuration (2), where the hyperbolic 
hedgehog is situated at the surface of the nematic drop. 
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FIG. 10. The free energy Fasa function of the position z\ 
of the water droplet for the configurations (1) and (2). For 
z\ > —3.5, (1) is stable, and (2) is metastable. The situation 
is reversed for —4.3 < zi < —3.5. Configuration (1) looses its 
metastability at z\ = —4.3. 

Moving the water droplet closer to the surface induces a 
repulsion due to the strong director deformations around 
the point defect. When the water droplet is placed far 
away from the south pole, i.e., at large zi, the dipole 
of configuration (1) forms and represents the absolute 
stable director field. At z\ = —3.5 the dipole becomes 
metastable but the system does not assume configuration 
(2) since the energy barrier, the system has to overcome 
by thermal activation, is much too high. By numerically 
calculating the free energy for different positions of the 
hedgehog we have, e.g., at z\ = —4.0 determined an en- 
ergy barrier of K 33 a » 1000 k^,T, where fee is the Boltz- 
mann constant, T the room temperature, and awl /im. 
At z\ = —4.3 the dipole even looses its metastability, 
the hyperbolic defect jumps to the surface at the south 
pole and the water droplet follows until it reaches its en- 
ergetically preferred position. On the other hand, if it 
were possible to move the water droplet away from the 
south pole, the hyperbolic hedgehog would stay at the 
surface, since configuration (2) is always metastable for 
Z\ > —3.5. The energy barrier for a transition to the 
dipole is again at least 1000 k-g,T. We have also inves- 
tigated the distance d\ of the defect from the surface of 
the water droplet. For z\ G [—2,4], d\ fluctuates be- 
tween 0.3 and 0.35. For z\ < —2, it increases up to 0.5 
at z\ — —4.3, where the dipole looses its metastability. 

D. Repulsion without Defect 

We return to the first geometry with two water droplets 
and homcotropic boundary conditions at all the surfaces. 
When we take either a uniform director field or randomly 
oriented directors as a starting configuration, our system 
always relaxes into the configuration sketched in Fig. O. 



FIG. 11. An alternative metastable configuration. Both 
droplets are surrounded by a —1/2 disclination ring which 
compensates the topological charge +1 of each droplet. An 
additional +1/2 disclination ring close to the surface of the 
nematic drop satisfies the total topological charge +1. 

Both water droplets are surrounded in their equatorial 
plane by a —1/2 disclination ring which compensates the 
point charge +1 each droplet carries with it |l^,^]. To 
obtain the total point charge +1 of the nematic drop 
there has to be an additional topological defect with a 
point charge +1. In the numerically relaxed director field 
we find a +1/2 disclination ring close to the outer surface. 
This configuration has a higher energy than the one with 
the hyperbolic hedgehog. It is only metastable. Since a 
transition to the stable configuration needs a complete 
rearrangement of the director field, the energy barrier is 
certainly larger than K 33 a sw 1000 k^T. We, therefore, 
expect the configuration of Fig. [ll] to be stable against 
thermal fluctuations. It would be interesting to search 
for it in the experiment. 

We use the configuration to demonstrate that even 
without the hyperbolic hedgehog the two water droplets 
experience some repulsion when they come close to each 
other. In Fig. |l2| we plot the free energy F versus the 
distance d of the two spheres. For large d the free en- 
ergy oscillates which we attribute to numerical artifacts. 
For decreasing d the free energy clearly increases, and 
the water droplets repel each other due to the strong de- 
formation of the director field lines connecting the two 
droplets. 
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APPENDIX: METHOD OF FINITE ELEMENTS 
FOR NEMATIC LIQUID CRYSTALS IN A 
CYLINDRICAL GEOMETRY 

The major advantage of the finite element method p^ ] 
is its ability to cope with arbitrarily complex geometries. 
In two dimensions, e.g., the integration area is subdivided 
into finite elements, which in the simplest case are trian- 
gles. In doing so, the boundaries of a complex geometry 
are well approximated by polygons. This is not possible 
within the method of finite differences, where the grid is 
defined by the coordinate lines. In our geometries strong 
director deformations occur close to the boundaries of 
the spheres. For a quantitative analysis it was, therefore, 
important to have these boundaries defined as smoothly 
as possible. 

To determine the tilt angle field the Oseen-Zocher- 
Frank free energy has to be minimized. The finite 
difference technique would directly address the Euler- 
Lagrange equation. On the other hand, the strategy 
of the finite element method is to start from a dis- 
cretized version of the free energy and then to apply 
a numerical minimization scheme. In our case, the re- 
duced free energy is written as a function F[0], where 
© = (0i, . . . , 6 n ) denotes the set of tilt angles on the n 
vertices of the net. It can be expressed as a sum over the 
free energies for the individual finite elements: 



three tilt angles at the vertices of the triangle. In order to 
derive Fj\Q^\ ®3 ] we start from the free energy 
of one triangle of area A(j): 



F i (6) = 27r / dzdppf[Q(p,z)} , (A2) 

JA( 3 ) 

The quantity /[0(p, zj\ is the Oseen-Zocher- Frank free 
energy density of Eq. (0), and an axial symmetry is as- 
sumed. Furthermore, in order to handle every finite ele- 
ment in the same manner, we introduce natural coordi- 
nates u, v: 



«W)(u, v) = zf + (z^> - z\ JI )u + (4 J; - z K l J >)v 



(A3) 



where pf^ and z\ JI are the cylindrical coordinates of the 
vertices of the j-th triangle (i — 1,2,3). In the new 
coordinates every triangle is right-angled and isosceles. 
The free energy now assumes the form 



Si) 



F S (Q) 



where 



2-jrAj x 

dv \ 
• < Jo 



dup u \u,v) f[Q^(u,v)} , (A4) 



A 



- f n U) 



Si)\ 



iP? 



(A5) 



is the Jacobian determinant of the coordinate transfor- 
mation. When the finite elements are small enough, the 
tilt angle field within one element is well approximated 
by a linear interpolation on the area of the triangle, 



Q u \u,v) 



9« 



(&i 3) -e[ j) )v 



(A6) 



We now insert (A6) and jAq ) into (A4). To simplify the 
integration over u and v we replace and p by their 
average values p$ = (p^ + p% + /4 )/3 and 8q = 

(6^ } + 9^ + 03 J) )/ 3 - This renders the integrals trivial, 
yielding the free energy of the j-th triangle in terms of 
the three tilt angles at its three vertices: 



W)7rflC?)i 



(A7) 



F[&] = J2F j [e[ j \ei j \e 



Uh 



a,U) qC?) a(i) 



G © 



(Al) 



The partial derivatives P and O z necessary to compute 
/[eW] from Eq. (§) follow from 



The sum runs over all I triangles, and Fj[Q^\ 6^] 

is the free energy of the j-th triangle as a function of the with 



(A8) 
(A9) 
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and 



,U) _ iJS> 



a) 



(A10) 



,<j) 



„(j) 
'l 

Jj) 



)/A„ 
p«)/A, . 



(All) 
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